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J> , Abstract 

o ■ 

The decay widths of the ++ and 1 + heavy- light hybrids to B(D) and pion are 
calculated by using the QCD sum rules. The interpolated current of the hybrid is cho- 
sen as gqj a G*T a h v (x). In order to simplify the calculation and avoid the ambiguity 
of three-point correlation function, a two-point correlation function between the pion 
and vacuum is used instead. The decay width of the ++ — > B(D) is about 12(16) 



(T) . MeV while the 1 + — > B{D) is around 0.4(1.8) MeV. We keep the leading order of 

\ \jM.Q expansion in our calculation for convenience. 

00 

1 Introduction 

It has been a long time for us to search the exotic hadrons such as the glueballs and 
hybrids. No candidate that has been found confirmedly so far. Recently, there are some 
evidences for the existence of hybrids resonance however. The E852 collaboration at BNL[|TJ 
has reported a J pc = 1~ + isovector resonance p(1405) in the reaction n~p — > rj^n with the 
mass 1370 ± I63J0 MeV and width 385 ± 40ii 05 MeV. Since the normal qq meson only has 
even spin in the channel {t]7i°), the state 1 h must be beyond the quark model. The Crystal 
Barrel collaboration has also claimed to find an evidence in pp annihilation which may be an 
resonance with a mass of 1400 ± 20 ± 20 MeV and a width of 310 ± 50±f MeVg]. Recently, 
the E852 collaboration put forth evidence for another J pc = 1~ + exotic meson p(1600)[[|, 
decaying to /wr, in the reaction ir~p — > 7t + 7t~tt~p, with a mass and width of 1593 ± 8 MeV 
and 168 ±20 MeV respectively. If these experiments have been confirmed in the future, they 
will provide a very strong evidence for the existence of the constituent gluon. 

Theoretically, the hybrids have been studied widely by various methods such as bag 
modelQ, flux-tube model[|J, QCD sum rules 0, lattice and some other models@. How- 
ever, there are few works about the heavy-light hybrids except for QCD sum rules 0] in 
which the spectrum of the heavy-light hybrid was given. Compared with bb and cc hy- 
brid, the heavy-light hybrid is easier to be dealt with since the heavy quark effective theory 
(HQET) may work in such a system. As we know, HQET has led to much progress in the 
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theoretical understanding of the properties of hadrons||. In such a framework, many phe- 
nomenology methods become easier to be controlled. For instance, QCD sum rule in HQET 
has been applied to various calculations of physical parameters, such as the decay constants, 
form factors and decay widths||. 

In the heavy-light hybrids case, the calculation in full QCD theory in Ref. |10| shows that 



the constituent of gluon gives a contribution more than l.OGeV to the mass, so the "light 
freedom" is too heavy and the availability of 1/Mq expansion seems problematic. Let's 
consider this problem. In the framework of 1/Mq expansion, the mass of heavy hadron can 
be expressed as 

M = M Q + m light _ freedom + \ 2 /M Q + (1) 

where the A have the dimension of mass and can be estimated by the mass split of the 
doublet such as (B ,B*) ,(D ,D*) because of the heavy quark spin symmetry. For instance, 
M B * — M B = 0.05 GeV, then A ~ 0.5GeV, which is roughly equal to mn ght _f reedom = 
Mb — M b = 0.5 GeV if we choose M b = 4.7GeV, therefore, the next-leading correction is 
around X/M b pa 1/10. Similarly, one can obtain A pa nought- freedom — M D — M c = 0.5GeV 
if one chooses M c = 1.3GeV, and A/M c pa 1/3 for doublet (D,D*). From the calculation 
in Ref. [|IIJ], the mass split of ++ and 1 h is 0.5GeV for b hybrid, so A = 1.6GeV ~ 



M b i-+ — M b and the next-leading correction is around X/m b pa 1/3. For the c hybrid, 



the mass split from Ref. [10] is O.SGeV, A = l.OGeV, which deviates a little more from 
mught- freedom = M c>1 -+ — m c pa 2.0GeV{it is also a little smaller than 1.6 GeV) , and the 
next-leading correction A/M c ~ 1.0. Comparing with B and D mesons , we can think it is 
still safe to use 1/Mq expansion in b hybrid system while it should be very careful to apply 
HQET to c hybrids. 



As to the decay property of hybrids, the results given by the QCD sum rules |TTJ gave 
a strong disfavor to the experiments. The decay width for some individual channel calcu- 
lated by them is much lower than the reported experimental results. As we know, all the 
predictions in sum rules before is dependent on the three-point correlation function, so the 
uncontrolled ambiguity resulted from the double Borel transformation will be brought in 
and the infrared divergence in the soft pion limit will appear too. To avoid these disadvan- 
tages, we make use of the two-point correlation function between the pion and vacuum in 
our derivation instead of the normal three-point correlation function. Besides, it simplifies 
our calculation. 

In this paper, we employ the HQET sum rules to calculate the decay width of the ++ 
and 1 _+ hybrids to B(D) meson and pion. For convenience, we keep only the leading order 
of 1/Mq expansion in our calculation, which may make some deviation in the c hybrid case. 
Our numerical results show that the decay width of ++ — > B(D) is about 12(16) MeV while 
the decay width of 1~ + -> B(D) is only 0.4(1.8) MeV. 

The paper is organized as follows. The analytic formalism of HQET sum rules for the 
decay of hybrid is given in Sec. 2. In Sec. 3, we calculate some pion's matrix element which 
is necessary in the sum rules for the calculation of decay width. The numerical results of 
the decay width were obtained in Sec. 4. We give the conclusion and discussion in the last 
section. 
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2 HQET sum rules for the decay of heavy-light hybrid 
mesons 



In this paper, we consider the following processes 

H b (0 ++ )(k) — > B(0-+)(k -q)+ n ± (q), (2) 
H b {r+){k) — > B(0-+)(k -q) + n ± (q), (3) 

where the Hb(0 ++ )(k) represents the hybrid with b quark and momentum k, which has 
jpc _ o ++ , Hb(l h )(&) represents the hybrid with b quark and momentum k, which has 
jpc _ i~+ Here, the electric charges of the mesons except for pion have not been written out 
explicitly. The cases of H c (0 ++ ) — > Dir^ 1 and H c (l~ + ) — > Dtc^ 1 are completely analogous. To 
calculate the decay widths of these processes, we consider the following two-point correlator 

A v (u',u,v)=i J dxe ikx (n ± (q)\T{j lv (x),j2(0)}\0) = A(u/,u)v v +B(u/,u)(-q u +q-w u ) (4) 

Where ji v (x) = gq^ ^G a ^ u T a h v {x) , ji{x) = h v j 5 q(x), h v is the heavy quark effective field with 
four velocity v, A(u',u) and B(uj',uj) are scalar functions of u and u', where uj and u' are 
defined as 

uj = 2k-v , u' = 2(k-q)-v. (5) 

A(uj',uj) and B(lj',lj) are determined through the spectral density saturated by the 
mesons corresponding to the interpolated currents, respectively. Here we mention that the 
interpolated current for the hybrid with a fixed J pc is not unique, different interpolated 
current corresponds to different state p. 

Before carrying out the operator product expansion, we can simplize the correlation 
function Eq. (|4|) firstly. Since the free heavy quark propagator in x representation in HQET 

is dr5(x — vt)^-^- ^and the interaction of the heavy quark with the gluon field A^ in the 
leading order of 1/Mq expansion is ghv ■ Ah. Then in the fixed-point gauge x^A^ = 0(which 
will be used throughout this paper), the full propagator of the heavy quark (0\T(h(x)h(0)\0) 
in the leading order of 1/Mq expansion is the same as the free one. Then the correlator (f|) 
reduces to the one including only light quarks. 
Performing the OPE in A u , we get 

A v = ij dxj k *{it ± {q)\T{j lv {x),32{U)M (6) 
= - J dxe i{k - q)x J rfr<5(x-t;r)rr{r(7r ± (g)|G^(0)g(-x)g(0)|0)} 
= -J dxe^Tririiir^lG^qqlO) - (tt^tv ■ DqqG^O)}}, 

where T = 7^ ^^75 and G^ w = igG a T a . In deriving the last equation in @, the quark field 
q{—x) is expanded around zero and only the first two terms are kept, the numerical results 
in section. 4 assure it a good expansion. We stress here that, since the perturbative theory 
does not break chiral symmetry d^j^ = 0, there is no perturbative contribution to the sum 
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rules in (|B|) in the chiral limit. The two matrix elements in (|G[) will be obtained through 



another suitable sum rules with the method used in Ref. |12|]. Their detailed expressions 
will be given in the next section. 

On the other side, A(ou',u) and B(uj',uj) can be represented with the spectral density 
through the dispersion relations. Because the perturbative contribution vanish, we suppose 
that the contribution of continuum states vanishes too. 

From the definition of Eq. (H), we have 

uj — uj' = 2q ■ v 

and one has to make the double Borel transformation on the two variables uj and uj' . The 
soft pion approximation was used in the spectral density in Ref. 0, which deduced func- 
tions of single variable. However, in the case of heavy-light hybrid's decay, the soft pion 
approximation is not good since the pion momentum may be large. Here we make a more 
reasonable approximation 

uj - uj' = 2q ■ v w 2(A - A') (7) 

in the infinite heavy quark mass limit, where A ~ m B — Mq and A' ~ m meson — Mq. Eq. 
(0) can be obtained from 

k 2 + q 2 - (k - qf = 2k ■ q » 2M Q q ■ v 

as Mq — > oo. We find that the difference between uj and uj' is around 1 — 2 GeV in the 
case of heavy-light hybrid. With the single pole terms and double pole term included in the 
spectral density, then A(u', uj) and B(lu', uj) are expressed as functions of the single variable 
uj', respectively 

At >\ — f H+ fm9H+mn'mj I+ m 2 n C {q\ 

A \ u ) - TWIT, . JvHT?, h 



{2h.'-uj') 2 M Q 2A'-u/' w 



(2A' - uj') 2 Mq- 2A 



uj 



where run, vn m and Mq are the mass of the hybrid, meson, and heavy quark, respectively. 
Co and c\ in the above equations are constants, fi are decay constants and gH ± m-K refers to 
decay amplitudes. For convenience, they are defined in full QCD theory as below 

(0\j H \H(0 + +)) = f H+ m 3 H+ K , (0|^|tf(l- + )> = f H -m%-e v , (10) 

(0\j D \D) = -if D m 2 D /M c , {0\j B \B) = -if B m 2 B /M b , 

(7r ± (g) J D|if(0 ++ )) = g H+Dn , (7T ± (q)D\H(l-+)) = g H - Dir e ■ q, 

(7r ± (q)B\H(0 ++ )) = g H+Bn , (n ± (q)B\H(l- + )) = g H - Bw e ■ q. 

Therefore, only one variable Borel transformation is needed. Taking the Borel transfor- 
mation 

~B?= lim (11) 

—uj' — > oo 
r = — uj' /t fixed 
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on both sides of the sum rules and differentiate the results in the variable 1/r, we get the 
expression for the decay amplitudes 

M 3 

9n +m , = q 4 2 [2A'A'(r) + A ]e^''\ (12) 

J H+ Jm"^H+" l m 
M 3 

9H-mn = , . Q 4 T [2Mff{r) + B ]e 2A ' /T , (13) 

where the A'(t) and B'(t) are the Borel transformed function of A(uj') and B(u'), respec- 
tively. Their expressions are obtained 

A! if) = 4V2{l3b 1 -2(m H -m m ) 2 d 2 ] + l(m H -m m )F 1 /4 + F 3 /l6]-}, (14) 

T 



B\t) = 4v / 2[2(m H -m m )rf 2 + F 1 /(12r)], (15) 

and 

A = 4v / 2[(m H -m m )F 1 /4 + F 3 /16], 
Bo - —, 

where the parameters and Fj will be given in the next section. 

Once the decay amplitudes are determined, it's straightforward to obtain the decay widths 
of the processes (EJ) and (H). 



+ (ei + e2^)g M g y ]75} a /3 ( - ) (16) 

ab 



3 The sum rules for pion's matrix elements 

The key to the sum rules is the determination of the two pion's matrix elements in ( 
The first pion's matrix element (vr ± (g)|G A1 j y gg|0) has already been given in Ref. 0(d) as 

(7^)1/^/^(0)^(0)10) = i{[g„u(ai + a 2 f) - iu^M + ie^Y q a ^b 2 

+ i{q^xu + quO\^)q x c 2 + i{q^o Xu - qucr\ fl )q x d 2 

'Ti 

where the definition and the numerical results of the coefficients can be found in the same 
Ref.. From the Lorentz covariance, the matrix element (7r(g) ± |(/} Q g)*gG j a 7 |0) can be written 

as 

(7T ± (g)| J D a ggG /37 |0) = {ifiq a cr^ - fzlakpl-y ~ Q-rlfi] + ihlaPfa ~ UiQapq-y ~ g^qp) (17) 

-highly ~ g a yl(i)}l5 + ij tfafrpQp ~ .Ma (<?/?7 7 ~ Q-ylp), 

where /j are some constants to be determined. By using the motion equation in the chiral 
limit and performing some special traces, we can obtain the following equations about f] 

h = -\Fi, h = U F i~ F *)> ^3 = ^3, (18) 

11 . , _ 1 



and ff vanishes. Where F\, F 2 and F 3 are some constants defined as 



(^(q^qGftvDaqlO) = ge tluap q p F 1 , (19) 
(7r ± (g)|g7 5 G Ml/ D a g|0) = ig(g au q p - g ap qu)F 2} (20) 
(^(^TsTmGVA^IO) = igF 3 . (21) 

In order to get the numerical results of the constants Fj, we take advantage of some 
other suitable sum rules, through which the constants Fj can be obtained by the following 
correlation functions, respectively 

K^M) = ^Jdxe^(0\T{qG^T a D a q(x),q l5 q(0)}\0), (22) 

n^(g) = tjdxe^iOlTiq^G^D^.q^qmiO), 

U%q) = zjdxe^^0\T{q l5lp G%T a D u q(x),qi P l 5 q(0)}\0). 

where the q or q refer to the u, d quark or antiquark corresponding to which in the tt ± . 

Keeping the matrix elements of OPE up to dimension six and omiting the radiative 
correction, we obtain the Borel transformed functions 

F 2 (M) = -F 1 (M), 

= {^(M 2 ) 2 n(0|?^|0)-^(0|«|0> 2 }fe— . 

where M is the Borel transform parameter of the sum rules. 

To obtain the numerical results, the constants and condensates are chosen as 

U = 132MeV, m n = UOMeV, m u + m d = 15MeV, (24) 

(0|— G 2 |0) = 0.012GW 4 , (0|gg|0) = -(0.24GeV0 3 , a s (0|gg|0> 2 = 8 ■ l(T 5 GeV 6 

7T 

and where the a s = 0.4. 

The numerical results are showed as Fig. 1 and Fig. 2, where one can find stable platforms 
in the region M ~ 1.1 — 1.6 GeV. they read 

F 1 = -F 2 = 4.0 • 10- 3 GeV 4 , F 3 = 2.3 • W~ 2 GeV 5 . (25) 



4 the numerical results of the decay width 

Before going on the numerical calculation of the decay width, there are still some param- 
eters necessary to be fixed. The heavy quark and meson masses are given as those in Ref. 
|T3| and the decay constants of B and D mesons are chosen as Ref. []14[ . The masses of the 



heavy-light hybrids have been computed in Ref.fllOf in the full QCD theory though the b 
quark mass(4.23 GeV) and some other parameters were chosen unsuitablely, we just recite 
the results of their calculation in this paper. 
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Table 1: Parameters input and the decay amplitudes. 





71 /T 

Mq 


77Z m 


p 
J 111 


m H {0 ++ ) 


m H (l +) 


b quark 


4.7 GeV 


5.28 GeV 


180 MeV 


6.8 GeV 


6.3 GeV 


c quark 


1.3 GeV 


1.87 GeV 


190 MeV 


4.0 GeV 


3.2 GeV 




/ff(0++) 




9H+m-K 


9H~mir 




b quark 


0.0012 


0.0012 


3.2 GeV 


1.2 




c quark 


0.004 


0.004 


2.0 GeV 


1.1 





The decay constants of hybrids hadn't been determined in Ref. [[Uj yet, so we fix them 
through the parameters and formulae given in it as 



M 2 ™R 

ft = ^e^IL(M 2 ), (26) 



where IL,(M 2 ) refer to the functions corresponding to the J pc = ++ and J pc = 1 + states 

3a r 1 r 1 m 2 
n s (M 2 ) = — |M 6 / dx / dyx(l - x)y 3 {3y - 2)e ^ (27) 

47T d JO JO 

! Ct rn 2 Ct ™^ 

— — (—G 2 )M 2 / dxx 2 e~ m^x -{mqq)M 2 / <ia;x(l — x)e~ m^ x 5 

8 7T J0 7T Jo 

IL(M 2 ) = ^|M 6 j dx C dyx{l - x )y 4 e~^ 

47T d JO JO 

I ct m 2 ct m 2 
+ —(—G 2 )M 2 dxx(2-x)e~u^ -(mqq)M 2 dxxil - x)e~^~* . 

24 7T JO 37T JO 

where wir is the mass of the hybrids, the Borel parameter M is chosen in the region where 
both the masses and the decay constants have platform and the scale of the ct s is set at the 
Borel variable. Their results are given in the table 1 too. 

With the help of Eq. (|i~2"|) and Eq. (0), the numerical results of the decay amplitudes 
are shown as Fig. 3 and Fig. 4. The platform is not good, we choose the value around the 
region r ~ 2.5 GeV, which is the suitable region in sum rules analysis, as our results. The 
numerical result shows that the next-leading contribution is less than 1/5 of the leading term 
at this region, so the OPE is well performed. All the parameters and some results calculated 
are collected in table. 1. 

As to the processes (Q) and (|3|) we considered, with the help of the formulae for two body 
decay, the decay widths in the leading order of 1/Mq expansion are given as 

r(lf(0++)->m(0-+)+7r)= = 9 \™ m2u ~f m i (28) 
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v v ; v ; ; 247T ^ 192vrm^ 

When all the numerical results are put in, the decay width corresponding to processes 
(§) and (H) are obtained 

T(H(0 ++ ) -> £((T + ) + tt) = 12Me^, T(H(0 ++ ) -> D(0~+) + tt) = 16MeV, (29) 
r(F(r + ) -> 5((T + ) + tt) = 0.4MeV, T(H(1~+) -> D(0~ + ) + tt) = 1.8MeV. (30) 

As we know, the decay of hybrids appears to follow the almost universal selection rule that 
gluonic excitations cannot transfer angular momentum to the final states as relative angular 
momentum, that is to say, the decay of hybrids to S + P final pairs of mesons are preferred, 
while the decay of hybrids to two S-wave mesons are suppressed [pl|. But the selection 



rule is not absolute, in the flux tube and constituent glue models, it can be broken by wave 
function and relativistic effects, and the bag model predict that it is also possible that the 
excited quark loses its angular momentum to orbital angular momentum ||16|| . Though the 
results obtained here are much larger than those calculated for the similar decay modes in 



the light hybrids case flllfl , which depended on the construction of three point correlation 
function, these results couldn't be excluded by theory. 

The reason that the decay widths of the 1 h hybrids are much smaller than those of the 
++ hybrids results from the additional factor ql/m 2 H in the width formulae (p8|) compared 
with that of the ++, s, since in the decay channel ifb jC (l~ + ) — > B(D)n, the final states are 
in a relative P-wave. 



5 Conclusion and Discussion 

In this paper, we discuss the decay of the heavy- light hybrids H(0 ++ ) and H(l~ + ) to 
B(D)ir by QCD sum rule approach in the framework of HEQT. With the help of the spectrum 
of the heavy-light hybrids flO||, we argue, in the b hybrid case, convergence of 1/Mq expansion 
is similar as that in the D meson case, while it is not so good for the c hybrid. The two-point 
correlation function instead of the normal three-point correlation function is employed for the 
derivation, and the heavy quark freedom is extracted out of the correlation functions in the x 
representation, which leads to the estimate of some pion's matrix elements and facilitates the 
calculation. In the leading order of 1/Mq expansion, the process of the calculation becomes 
much easier. 

To find the hybrids, one approach is to look for an excess of observed state over the 
number predicted by the quark model, the other is to search for quantum numbers which 
cannot be accommodated in the quark model, so the experiments place main efforts on 
the exotic mesons' searching and the current experiments concern mainly the light hybrids 
sector, however, it's necessary to extend the search to the heavy-light hybrids sector. The 
calculation shows that the decay width of H(0 ++ ) — > B(D)ir + is around 12(16) MeV, while 
the decay width #(1~ + ) -> B(D)tt + is only 0.4(1.8) MeV. If people think the total width 
of hybrid is around 200 MeV, the channel B(D)ir may not be dominant in hybrid's decay, 
for instance, the channels B, D(B* , D*)rj' may be larger because of QCD anomaly. The 
calculation of these and other channels is beyond this paper. 
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There are some uncertainties unavoidable in this paper, which result from such as pole 
approximation, 1/Mq expansion(we keep only the leading order term) and OPE(we keep only 
the first two terms of OPE in (|6|)), omitting the higher orders would give some uncertainties. 
Besides, the calculation is dependent on the numerical results of the masses and decay 
constants of the hybrids taken from the Ref . fl(| , which were determined by some unpopular 



parameters. To follow the calculation of Ref. [10| with the new values of the parameters is 



beyond the goal of this paper. Certainly, it's necessary to identify the spectrum and decay 
width of heavy-light hybrids through other models and other sum rules, for example, light 
cone sum rules. Only when most of the models give detailed research on hybrids, can we 
declaim that we have a good understanding about the hybrids. 
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Figure caption 

Figure 1: F\ versus Borel variable M. 

Figure 2: F3 versus Borel variable M. 

Figure 3: Decay amplitude of ++ heavy-light hybrids versus Borel variable r. 

Figure 4: Decay amplitude of 1~ + heavy-light hybrids versus Borel variable r. 
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